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0-{\ ABSTRACT 

Q . We present the power spectrum of the fluctuations in the cosmic microwave back- 

' ground detected by the Very SmaU Array (VSA) in its first season of observations in 

I its compact configuration. We find clear detections of first and second acoustic peaks 

. a.t £fii 200 and £^ 550, plus detection of power on scales up to ^ = 800. The VSA 

power spectrum is in very good agreement with the results of the BOOMERANG, 
DASI and MAXIMA telescopes despite the differing potential systematic errors. 
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1 INTRODUCTION 

Anisotropies in the Cosmic Microwave Background (CMB) 
have now been detected b y many ex perimen ts (most recently 
iNettcrficld ct al.' ( 200l\ Le e et al.l J200ll) . iHalverson et all 
(2002), Padin et al. (2001)). At present the most successful 
model for explaining the origin of these fluctuations postu- 
lates that they are seeded in the very early universe by quan- 
tum perturbations with random phase, followed by a period 
of inflationary expansion. The fluctuations in the CMB are 
predicted to be Gaussian in nature, and hence can be com- 
pletely characterised through their power spectrum. A fur- 
ther prediction is that the power spectrum will show acoustic 
peaks due to plasma oscillations on scales smaller than the 
sound horizon at the surface of last scattering. 

In this paper we present the power spectrum of the 
CMB fluctuations detected by the Very SmaU Array (VSA) 
on spherical harmonic modes £ «150-900. We outline how 
the fully calibrated time-stream data are converted into a 
power spectrum and the various data checks that we have 
performed to confirm the validity of our analysis. This pa- 



per is the third in a series of four papers which report the 
results of the flrst season of observations m ade using the 
VSA in its compact configuration. Paper I iWatson et alJ 
l2002h describes the design of the VSA and our experimen- 
tal method; the observational strategy, foreground removal 
and reduction techniques for t he data analysed in this pa- 
per are described in Paper II iTavlor et aljl20o3) : finally, 
the cosmological implication s of the VSA power spectrum 
are discussed in ■Rubifio-Martin et alJ 120021 (Paper IV). 



2 OBSERVATIONS AND INITIAL DATA 
PROCESSING 

2.1 The observations 

The Very Small Array (VSA) is a 14-element interferometer 
array for Cosmic Microwave Background (CMB) observa- 
tions; it operates at a frequency between 26 and 36 GHz with 
a receiver bandwidth of 1.5 GHz. In its compact configura- 
tion, used here, the instrument is most sensitive to angular 
structure in the range £ «150-900. As well as the 14-element 
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array, there is also a single baseline interferometer used for 
radio source flux measurements. A fuller description of the 
instrument is given in Paper I. 

The present observations were made at a frequency of 
34 GHz in the period 2000 September to 2001 September 
and were centred on three separate areas of sky. Each VSA 
observation maps a region with FWHM of 4.6°. The pro- 
cess of field selection is discussed in Paper II; selection was 
based primarily on low Galactic (synchrotron, free-free and 
dust) emission and an absence of known bright foreground 
sources. Overlapping fields were observed in each area in or- 
der to reduce the sample variance, to increase the resolution 
in ^-space and to allow for direct assessment of data relia- 
bility and the detection of any residual instrumental effects 
(Section E^. The array configuration used was designed to 
provide an approximately uniform spread of interferometer 
baselines while retaining a reasonable aperture filling factor. 

The amplitude and phase calibrations of the individual 
interferometer baselines were checked both by frequent short 
meeisurements of a rmmber of secondary calibration sources 
and also by regular longer observations of a small number 
of primary calibration sources. The overall calibration pro- 
cedure is described fully in Paper I; the overall accuracy of 
the calibration in flux density and temperature is 3.5% . 

An important aspect of the VSA is the inclusion of a 
separate interferometer, comprising two 3.7-m dishes on a 
9-m north-south baseline, for determining the flux densi- 
ties of foreground sources (radiogalaxies and quasars) falling 
within the observed fields. The positions of all sources which 
might affect the VSA observations were obt ained from sur- 
vey o bservations using the Ryle Telescope dWaldram et alJ 
120021) : these positions were subsequently observed concur- 
rently with the main VSA observations in a series of reg- 
ular pointed observations and the measured fiux densities 
subtracted from the visibilities observed by the main array 
(Paper II). 

2.2 Initial Data Processing 

The data are calibrated and processed as described in Papers 
I and II. Early tests of the telescope revealed an unwanted lo- 
cal spurious signal, particularly evident on the shorter base- 
lines. The procedures for removing this signal from the data 
are described fully in Paper II; the tests presented in Sec- 
tion 14.21 below indicate that the removal of this signal is 
complete. Similar filtering processes were used to eliminate 
signals arising from the Sun and the Moon. 



3 DERIVATION OF THE POWER SPECTRUM 

Since the three sets of mosaiced fields are widely sepa- 
rated from one another on the sky, the data correspond- 
ing to each field are analysed individually. The deriva- 
tion of the GMB power spectrum for each mosaic is per- 
formed using the maximum -likelihood method presented in 
iHobson fc Maisingerl (|20o3), which we summarise here. The 
results from each mosaic are then combined, as outlined be- 
low, to produce the final estimate of the power spectrum. 

Since the number of individual visibility measurements 
for each mosaic is very large (~ 800,000), it is first neces- 
sary to compress these data in some way. For each separate 



field in the mosaic, the visibilities are binned into cells in the 
uw-plane; for the case in which the instrumental noise covari- 
ance matrix is diagonal, this corresponds to the maximum- 
likelihood solution for the value of the binned signal visibil- 
ity in each cell. This is analogous to the 'map-making' step 
in the analysis of single-dish GMB observations, in which 
tim e-ordered dat a are binned into pixels on the sky (see, 
e.g. lBorrilll ^1999^ ) . Since we are not interested here in mak- 
ing accurate GMB maps from the binned visibility data, 
the uw-plane is simply divided into equal-area cells of size 
Alt = 3 wavelengths. As the aperture function of the com- 
pact VSA is well-modelled by a Gaussian with a FWHM of 
12 wavelengths, this cell size ensures that the uw-plane is 
comfortably over-sampled, while reducing the total number 
of visibility measurements for each field significantly (to ~ 
2500). 

The binned visibilities are thus the basic input to the 
likelihood analysis for the GMB power spectrum. The com- 
pact VSA is sensitive to the multipole range I ~ 150 — 900, 
and the effective aperture function (after mosaicing) has a 
FWHM of M = 83. We thus divide the total ^-range into 10 
spectral bands each of width A^, in order to limit the cor- 
relation between adjacent bands. In each band, we assume 
l{i + l)Ci to have a constant value Cb {B = 1,2,..., 10). 
These flat power levels in each band are then the parameters 
to be determined in our likelihood analysis of the data for 
each mosaic. We denote these band powers collectively by 
the parameter vector a. 

Assuming the GMB emission and the instrumental noise 
to be Gaussian random fields, the log-likelihood of obtaining 
the binned visibility data vector v (which contains the real 
and imaginary parts separately of each binned visibility), 
given some set of flat band powers a, is given by 

ln£ = constant — i [ln|C(a)| +w'C~^(a)w] , 

where C{a) = S{a) -I- iV is the sum of the predicted sig- 
nal covariance matrix and the noise covariance matrix. The 
maximum-likelihood GMB power spectrum a is calculated 
using a simple iterative numerical maximisation algorithm. 
Starting from an initial guess ao (which is unimportant), in- 
dependent line maximisations are performed for each band 
power Cb in turn, while keeping the others fixed. The whole 
solution vector a is then updated and the process repeated 
until convergence is obtained. This typically requires around 
5 iterations. 

The well-defined correlation structure of visibility data 
in the iiu-plane allows each line-maximisation to be per- 
formed using only the subset of visibilities that are sensi- 
tive to the band power being varied, thereby speeding-up 
the evaluation of the likelihood function, while keeping the 
calculation exact. For a single VSA mosaic, the maximum- 
likelihood solution can be obtained in around one hour on 
an 8-node Beowulf Cluster with 1.8-GHz AMD Athlon pro- 
cessors. 

The uncertainties in the derived maximum-likelihood 
GMB power spectrum are estimated in two complementary 
ways. Assuming the likelihood function in the parameter 
space a to be well-approximated by a multivariate Gaus- 
sian near its peak, the covariance matrix of the parameter 
uncertainties is given simply by (minus) the inverse of the 
curvature (or Hessian) matrix at the peak d. This matrix 
is easily evaluated numerically in a few hours of GPU time. 
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The square-root of the diagonal elements of the resulting co- 
variance matrix give the standard error ACb on each band 
power, whereas the off-diagonal elements provide a measure 
of the correlation between the band power estimates in dif- 
ferent spectral bins. We find, typically, that adjacent bins 
are anti-correlated at around the 5 — 20% level, and for more 
widely separated bins the correlation is negligible. 

Using the covariance matrix clearly produces symmet- 
ric error bars on each band power Cb, which may be a poor 
representation of the uncertainty, especially for poorly con- 
strained band powers. An alternative approach is to make 
use of the fact that the band power estimates in different 
spectral bins are only weakly correlated, so that the off- 
diagonal elements of the covariance matrix are small com- 
pared to those on the diagonal. In this better rep- 
resentation of the uncertainty in each estimate may be ob- 
tained by directly evaluating the likelihood function through 
the peak d, along each parameter direction in turn. In the 
ideal case, where the band power estimates are independent, 
the resulting curves would be the marginal distributions of 
each band power. Since the likelihood function can be eval- 
uated very quickly along each direction in parameter space, 
the resulting 'marginal' distributions can be calculated for 
each VSA mosaic in around 2 hours, and provide a useful 
cross-check of the standard errors obtained from the covari- 
ance matrix at the peak. We find that for spectral bins in 
which the flat band power is tightly-constrained, the shape 
of the likelihood function is very close to Gaussian. How- 
ever, for bins in which the power level less well constrained, 
the likelihood function is better described by an offset log- 
normal distribution; this is discussed further in Paper IV. 

Once the likelihood functions for the flat band power 
Cb in each spectral bin have been obtained separately for 
each VSA mosaic, the mosaics are combined simply by mul- 
tiplying together the respective likelihood functions in each 
bin. This assumes that each mosaic provides an independent 
measurement of the CMB power spectrum in each spectral 
bin, which is valid given that the three sets of mosaiced 
fields are widely separated on the sky. The correlation be- 
tween the resulting band power estimates in different bins 
is easily obtained by calculating the covariance matrix, as 
described above, at the new joint optimum d. 

All the above functio ns are implemented using the 
Madcow analysis package (iHobson fc Maisingei^l2002^ . 

3.1 Windovir functions 

Because of various instrumental effects (e.g. non-uniform uv 
coverage and the finite size of the primary beam of the tele- 
scope), a bin samples the underlying power spectrum C; 
through a window function W{1). For the _Bth bin, the mea- 
sured power corresponds to 

PB = 2^ — 

I 

where Ce = £{£ + l)Cf/(27r). 

The window function for a given bin is determined as 
follows. The range of Fourier modes to which a given vis- 
ibility is sensitive are given by the Fourier transform of 
the primary beam for a s ingle field observation (see e.g. 
iHobson fc Maisingeil (l2002l) 'l. For a Gaussian primary beam. 
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Figure 1. Window functions for the combined data set. The func- 
tions are normalised to unit area, and different bins are plotted 
with different linestyles to allow easier visual differentiation. 



^(u) = 27r(T^ exp(— 27r^(T^|u|^), 

where u is the vector in the visibility plane (measured in 
wavelengths). For mosaiced observations, the effective beam 
is a superposition of displaced primary beams. We can think 
of this superposition as a convolution of a centred primary 
beam with a sum of delta functions at the beam centres. In 
the Fourier domain this corresponds to a multiplication of 
the aperture function with the 'Fraunhofer diffraction pat- 
tern' of the beam centres 

Agfj(u) = E,A(u)c^---^ 

where Xj are the positions of the beam centres from the ori- 
gin, measured in radians. Note that in general this aperture 
function is complex. 

For a given bin, the weighted complex sensitivity map 

is 

•S'(u) = WkA^f^iu - Vfe), 

k 

where index k runs over visibilities in a given bin. These 
visibilities have instrumental weights Wk and positions on 
the visibility plane v^. 

The un-normalised window function is then given by: 

/ — 2-77 
|S(^,0)|'d0, 

where £ — 27r|u[. The normalisation constant A/" can be triv- 
ially found using J^^{2£ + l)W{£)/{2£{£+ 1)) = 1. The win- 
dow functions for the combined VSA data set are plotted in 
Figure Q showing the small degree of correlation between 
adjacent bins. These functions are used in Paper IV in the 
estimation of cosmological parameters. 



4 DATA CHECKS 

The complete process of editing and filtering the data and 
also the subsequent stages of data reduction were carried 
out independently by the Cambridge group and, jointly, by 
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the combined lAC and JBO teams. A comparison of the two 
sets of results showed good agreement, the effect of any dif- 
ferences being small compared to the intrinsic uncertainties 
on the final power spectra. 

4.1 Test of data reduction procedure 

Aspects of our data reduction procedure, such as filtering 
and calibration, could potentially have introduced system- 
atic errors into the VSA data. In order to test this, we pro- 
duced a realisation of the CMB sky and simulated a mock 
VSA observation including such instrumental effects as vis- 
ibility quadrature errors, phase steps due to path compen- 
sation and thermal noise. We analysed these simulated data 
using our standard reduction procedure and then produced 
a power spectrum as described in section |3 We found this 
to be entirely consistent with our input model and so con- 
cluded that our method of data reduction did not introduce 
any significant systematic errors. In particular, it demon- 
strates that any correlations in the visibilities introduced by 
the filtering process has not observable effect on the derived 
power spectrum. 

4.2 The Local Spurious Signal 

The local spurious signal, fully discussed in Paper I, was 
found to depend only on the antenna tracking angle, and not 
on the table elevation. Therefore, we would expect the spu- 
rious signal to be identical for different fields with the same 
declination observed over the same hour angle range. This is 
easily confirmed by combining the (unfiltered) baseline time- 
series of two such fields by both addition and subtraction. 
We find that adding the two fields enhances the spurious 
signal, whilst subtraction entirely removes it. Whilst this 
technique of addition and subtraction is adequate for de- 
tecting the presence of spurious signal in unfiltered data, it 
is insufficiently sensitive to detect possible low level residual 
signal once filtering has been applied. 

The increase in detection sensitivity that we require 
in order to test the filtered data for residual spurious sig- 
nal can, however, be obtained using a modified version of 
a MEM al gorithm used for extra ction of CMB signal from 
VSA data jMaisinger et al.lll997r) . 

We add an extra term to the MEM reconstruction which 
is the signal that is identical in the CMB datasets. We then 
consider the case of two fields at identical declination. As 
the CMB signal in the two fields will not be identical, and 
the noise is random in each case separately, any common 
component to the two fields will be spurious signal. 

We tested this algorithm by applying it to pairs of simu- 
lated observations of CMB fields assuming CDM primordial 
fluctuations, to which we added an identical component with 
an rms level such that it was not the dominant signal. We 
used a variety of common components, including a scaled 
down version of the unfiltered spurious signal. We found 
that the MEM algorithm was able to reconstruct the com- 
mon signals well, recovering the structure excellently, and 
recovering the amplitude of the signal to within 10%. As 
we are primarily interested in whether or not the spurious 
signal is still present, as opposed to any accurate quantifica- 
tion, this is perfectly adequate. Note that, even for entirely 
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Figure 2. Power spectra of the recovered common signal from the 
two VSA2 fields (squares) and from a pair of simulated datasets 
with no common signal (i.e. independent noise only) (crosses). 
The pairs of points have been separated laterally for clarity. Cov- 
erage of the uv-p\ane and the thermal noise level are identical for 
both the real data and the simulation. 

independent fields, the reconstructed shared signal power 
spectra are not zero (Figure I^J, but rather show a value in- 
creasing as consistent with the correlation between two 
white noise signals. 

The two VSA2 fields had identical declination to allow 
testing of our filtering procedures. The MEM algorithm was 
applied to these fields, and the results compared with CDM 
realisations with identical uv-coverage and thermal noise in 
which no shared signal is present. The results (Figure |5J 
show that the common signal found from the real data fields 
and the simulations agree within the errors, the difference 
between the two sets of points being less than 4% of the 
measured CMB for values of ^ < 750. In the same manner, 
we find no evidence for residual spurious signal in the pairs 
of VSA1,3 fields at similar declinations. 

Non-Gaussianity analysis of the binned visibilities al- 
lows us to locate and remove the few remaining visibilities 
contaminated with spurious signals down to a low level. The 
removal of these points has a negligible effect on the final 
power spectrum, giving us confidence that we are subtract- 
ing the spurious signal to a level well below that which could 
affect our results. The full details of this analysis will be 
published in Savage et al. in prep. 

4.3 Foreground Source Subtraction 

Radio galaxies and quasars are a significant contaminant 
of the CMB at microwave frequencies and in the higher i- 
ranges will dominate the CMB signal, making it essential to 
remove their contribution. Tests with simulated fields have 
been used to assess the potential contribution, before sub- 
traction, of these sources to the final CMB spectrum. We 
generated ten realisations of the CMB sky using a partic- 
ular CDM model and added to these the sources that we 
have observed to be present in the VSAl field. We have then 
compared the power spectra recovered from simulated VSA 
observations with and without sources. Two of these simu- 
lations, representing the range of results obtained from the 
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Figure 3. Two examples of recovered power spectra from simu- 
lated CMB observations with (filled circles, solid lines) and with- 
out (open circles, dotted lines) the known sources in the VSAl 
field added. The differences arise from chance interactions be- 
tween the sources and individual CMB features. 

ten simulations, are shown in Figure |21 It is apparent that, 
with the fairly small number of sources in any one VSA field 
(typically 12 sources), the impact on the power spectrum is 
unpredictable. Although the main contribution tends to be 
in the higher ^-bins (with errors potentially reaching ±100% 
at ^ = 900), changes of up to ±10% can occur in the lowest 
^-bin; these are due to the chance superposition of interfer- 
ence fringes. 

To determine the effect of residual sources on the VSA 
results, two further simulations without a CMB contribution 
have been carried out. We base these simul ations on the 15- 
GHz source counts from lTavlor et all i200ll) and extrapolate 
up to 34 GHz using a mean spectral index of a = 0.55. The 
first (Figure|l]upper plot) comprises the contributions of the 
known point sources and a statistical distribution of weaker 
sources; the second (lower plot) includes only the statisti- 
cal contribution of weaker sources and gives an indication of 
the possible residual contribution to the VSA power spec- 
trum after subtraction of known sources. It is clear that, 
with no source subtraction, the CMB data are significantly 
compromised for £- values ^ 600. After subtraction of known 
sources, the contribution is reasonably small for ^-values up 
to about 1000. As demonstrated by the results of Figure |21 
the contribution of sources to the observed spectra can not 
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Figure 4. Simulated power spectra for known sources plus a 
statistical distribution of weaker sources (upper plot) and for the 
distribution of weaker sources alone (lower plot). For comparison 
the dashed curve shows the predicted power spectrum for a CDM 
model. 

generally be predicted from a simple combination of sim- 
ulated power spectra, although such an approach becomes 
feasible in the limiting case of many weak sources per syn- 
thesised beam. 

We can also estimate the residual source contribution 
to the power spectrum using the preliminary 34 GHz source 
count derived from the source subtractor observations in Pa- 
per II. Integrating the count from zero flux to our complete 
source subtraction limit of 80 mjy and converting to units of 
AT /To, the source power is given by Ce = 7.7 x 10~^®, corre- 
sponding to a power spectrum value of ^^"^^Ce = 580 fiK^ 
at £ = 800. For the case with effectively no source sub- 
traction, taking the upper flux limit to be 0.5 Jy, we find 
Ce = 3.7 X lO"^'^, equivalent to 2800 ^iK^ at ^ = 1000; these 
results are in good agreement with the extrapolations from 
the 15 GHz counts. 

The precise residual contribution due to sources is sub- 
ject to uncertainties in our knowledge of the true weak- 
source distribution; the values shown in Figure 2] are likely 
to be over-estimates since some sources with fiux densities 
less than 80 mJy were actually subtracted. In order to assess 
the possible impact of these residual sources on derived pa- 
rameters (Paper 4) the cosmic parameter analysis has been 
repeated after subtracting the residual contributions shown 



6 Paul F. Scott et al. 



Figure 2] from the observed power spectrum. The changes 
in the fitted parameters are small, the largest effect being a 
reduction of 0.05 in ris, representing a change of about 0.5 
s.d. 



4.4 Galactic Foregrounds 

The diffuse Galactic foregrounds are discussed in Paper II, 
where it is suggested that the total contribution from the 
three components of the Galactic foreground (synchrotron, 
free-free and possibly spinning dust) amount to no more 
than about AT* = 5 /iK at an angular scale of 1° {I ~ 200). 
Here we discuss their effect on the results presented in 
this paper. The free-free and synchrotron components are 
relatively well-known and are each expected to contribute 
about AT = 1 — 2 /^K. The spinning dust component is 
more uncertain, but may be the dominant component, per- 
haps contributing up to AT ^ 5 /iK. This is based on a 
dust-correlated component with a correlation coefficient of 
10 MK/(MJy sr"^) at 100 /im. This coefiicient is stiU to be 
clearly demonstrated observationally. The power spectra of 
the diffuse foregrounds falls with increasing I (for example, 
[Giardin o et al. (2001)) and is likely to be half these values 
at £ - 500 . 

Any Galactic contribution adds in quadrature with the 
CMB signal and hence at the position of the first CMB peak 
[l ~ 200) which has ATcmb ~ 75 /iK, 5 /iK of foreground 
signal will increase the observed signal by 0.17 /^K. Simi- 
larly, at f ~ 500 where ATcmb ~ 45 /kK, the increase will 
be 0.28 /iK at the most. We see therefore, for the VSA fields, 
and at a frequency of 34 GHz, the contibution from Galac- 
tic foregrounds is likely to be negligible (< 1 percent). A 
similar conclusion, albeit based on observations of different 
regions of sky, was reached by Halverson et al. (2002) who 
found that Galactic emission made a negligible contribu- 
tion to the observed CMB power spectrum. A more com- 
plete cross-correlation analysis to investigate the contribu- 
tion from dust-correlated emission is in progress (Dickinson 
et al. in prep.). 

4.5 Noise Estimation 

The likelihood analysis used for power-spectrum estimation 
requires an accurate estimate of the rms noise level. Since 
on individual baselines the contribution of the CMB to the 
individual data samples is very small, the noise level can be 
obtained directly from the standard deviation of the data, af- 
ter the filtering and fiagging processes have been completed. 
The noise level associated with each of the binned visibili- 
ties (Section 3) is obtained from an appropriately weighted 
combination of the noise levels of the data in each bin. As an 
additional check, the scatter in the data points contributing 
to each bin has also been used to provide a noise estimate. 
Consistent estimates of the noise level were obtained by the 
two methods. The overall noise estimate is accurate to 2.5%. 

The sensitivity of the likelihood analysis to errors in the 
noise level has been tested by analysing the same dataset 
with different assumed levels of noise; the corresponding 
changes in estimated CMB power produced by the above 
uncertainty in estimated noise level are less than 1% for i- 
values up to ~ 750. 
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Significance 


VSAl 


977 


1033.7 


0.10 


VSAIA 


1884 


1947.7 


0.15 


VSAIB 


1387 


1410.3 


0.33 


VSA2 


915 


984.2 


0.06 


VSA2-OFF 


1384 


1420.2 


0.24 


VSA3 


1287 


1356.2 


0.09 


VSA3A 


2003 


2094.2 


0.08 


VSA3B 


1584 


1660.1 


0.09 



Table 1. The values for data splits on each of the VSA fields. 
In each case the visibility data from each field was split in two 
according to epoch and the of the difference vectors formed. 
Also tabulated are the number of degrees of freedom (DOF) and 
the significance of each value; the significance is given as the 
probability to exceed the observed value in the x^ cumulative 
distribution function. 

4.6 Tests 

As an additional check on the consistency of the data, we 
computed the statistic for a variety of splits on the 
visibility data. The data from each VSA field were split 
into two and the difference vector, A = {v\ — V2), formed 
where vi and 112 are the visibility vectors for each half 
of the data. The statistic was then calculated using 

— A*Cj^^A, where Cn is the noise covariance matrix 
given by Cn = (Cni -f Cn2). Since for each field the sky 
signal measured in each of the two halves of the data should 
be the same, the statistic can be used to test for the pres- 
ence of systematic errors in the data. The data for each VSA 
field were split in two according to observing epoch and the 

values and associated significances are given in Tabl(0 
The data split on the VSA3 field also corresponds to a split 
between day and night observations. This was the only field 
to be observed in this way, and its x^-value confirms the 
consistency of daytime versus nightime observation. 

The consistency of the power spectra derived from each 
of the 3 VSA mosaiced fields (VSA1M,VSA2M and VSA3M) 
was also compared by forming the x^ statistic on pairs of 
power spectra. In this case the x^ value is given by 

= (ai - a2f{H^^ -f H2^y^{ai - 02), 

where oi and 02 are the two sets of bandpowers and H^^ 
and if are the corresponding inverse Hessian matrices. 
Since x^ statistic assumes that the likelihoods are Gaus- 
sian, we use the Hessians calculated using an offset log- 
normal approximation (see Paper IV for further details). 
The x^ values (and significances) for the VSA1M/VSA2M, 
VSA1M/VSA3M and VSA2M/VSA3M power spectra com- 
parisons are 7.6 (0.67), 9.82 (0.46) and 5.03 (0.89) respec- 
tively. In each case there are ten degrees of freedom in the 
power spectrum analysis. 



5 THE POWER SPECTRUM 

The filtered and source-subtracted data for each of the VSA 
field s have been analysed using t he Madcow software pack- 
age iHobson fc MaisingeJ 1200a) as described in Section 3. 
The combination of the three mosaiced spectra is shown 
in Figure |^ The bin width used is /S.I — 83, which gives 
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Table 2. The power spectrum from combining the three VSA 
fields. The two sets of bin numbers (1, lA etc) refer to the main 
and alternate binnings, the latter being shifted by half a bin 
width. All the bins have A£ = 83. The reported errorbars corre- 
spond to 68% confidence limits and were calculated by enclosing 
68% area under a likelihood curve assuming independent bins. For 
further data analysis it is neccessary to use full window functions 
and covariance matrices that can be downloaded from our web- 
site. In addition to these errors, there is an overall 7% calibration 
uncertainty in power. 



B 


Cb,B-2 


Cb,b-i 


Cb.b 


Cb,b+i 


Cb,B+2 


Covs^s 


1 






1.00 


-0.06 


0.00 


18.78 


2 




-0.06 


1.00 


-0.13 


0.00 


15.53 


3 


0.00 


-0.13 


1.00 


-0.07 


0.02 


4.86 


4 


0.00 


-0.07 


1.00 


-0.18 


0.06 


1.98 


5 


0.02 


-0.18 


1.00 


-0.18 


0.01 


2.46 


6 


0.06 


-0.18 


1.00 


-0.15 


0.05 


6.19 


7 


0.01 


-0.15 


1.00 


-0.22 


0.02 


5.96 


8 


0.05 


-0.22 


1.00 


-0.16 


0.03 


10.06 


9 


0.02 


-0.16 


1.00 


-0.21 




12.41 


10 


0.03 


-0.21 


1.00 






35.88 



sensibly sample the power spectrum of our data. Numerical 
values for both binnings are given in Table |5| 

The plotted error bars contain the contributions from 
both thermal noise and sample variance, but not calibration 
errors, which introduce a completely correlated uncertainty 
in all the points of ±7 percent. Errors from pointing and 
primary beam uncertainties are negligible. Since the tem- 
perature sensitivity of the VSA compact configuration falls 
off dramatically after i = 800, all data above this have been 
binned together. 



6 DISCUSSION 

6.1 The VSA power spectrum 

The power spectrum shown in Figure |5] shows a clear de- 
tection of the first peak at £ ~ 220, and power at the level 
of about 2000 ^iK^ between £ = 300 and £ = 900. We have 
attempted to quantify the detection of a second peak at 
£ ~ 550, as this is the region of the power spectrum with the 
largest anti-correlations between adjacent bins (see Table|^; 
bin 6 centred at £ — 556 is anti-correlated with its neigh- 
bours at the ~ 20 percent level. We made Monte Carlo sim- 
ulations of the five Ce points between £ = 390 and £ — 722. 
Points were draw from the 5-dimensional Gaussian distri- 
bution described by the actual correlation matrix of these 
points, but with mean values e qual to the weighted mean o f 
the five actual points. Following lHobson fc Maeueiiol l)l99fll . 
who define a normalised convexity about a power spectrum 
point, we calculated the change in normalised convexity of 
the power spectrum about the inner 3 points, defined by 




where C = {C4 + Cs-3{C5 + C7))/2 + 2C6 and a is the overall 
error in C given from the errors in the individual points ai by 
o-^ = {al + al + 9{al + a^))/4+4a'i. This effectively compares 
the hypothesis that the power spectrum in this region is fiat, 
to the one that it is described by 'trough-peak-trough'. 

In 1000 realisations, we found only 27 instances of E 
being larger than the value observed in the real data of 
E = 2.2. We therefore conclude that the observed second 
peak is detected at 97 percent confidence. The power spec- 
trum is completely consistent with the adiabatic infiationary 
models, fits to which are discussed in Paper IV. 



Table 3. The correlation matrix dj for the combined VSA data 
set (main binning only) . These are calculated by normalising the 
covariance matrix of errors. Note that the correlation is only sig- 
nificant for adjacent bins. The values of matrix for which 
is not reported can be assumed to be zero. The final column 
gives the diagonal elements of the covariance matrix in units of 
10^ X /iK* 

weakly-correlated errors in each bin. The actual correlations 
between bins are given by the correlation matrix (Table I^J- 
To reduce the bias in assessing features in the power spec- 
trum caused by the settings of the bin centres, we have also 
calculated the power spectrum with bin centres shifted by 
one half a bin width to the right of the original bin centres. 
These results are shown in Figure|^with dashed error bars. 
Adjacent 'double- binned' points are highly correlated but do 



6.2 Comparison with other experiments 

In Figure |S| we compare the new VSA power spectrum 
plotted with those from BOOMER ANG jNetterfield et al. 
2001), DASI jHalverson et al.ll2002l) . MAXIMA jLee et al, 
2001D . Only single-binned (weakly-correlated) points are 
shown. We have attempted to compare the random and cor- 
related errors on the various experiments in a consistent 
way, difficult though this is on a single plot. Two sets of 
error bars are shown for each plot; the smaller bars indicate 
68 percent confidence limits from the random (thermal and 
sample variance) errors, while the larger error bars repre- 
sent systematic (calibration and beam) errors as reported - 
all the points from a single experiment are able to move up 
or down within the larger error bars. It is significant that, 
although the overall errors of the different experiments are 
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comparable, the relative contribution of systematic errors 
(including beam uncertainty) is much smaller for the inter- 
ferometric data. 

The agreement between the experiments on the exis- 
tence, heights and positions of two peaks and of power at 
higher I, is evident. This is particularly significant given the 
very different experimental techniques involved and the dif- 
ferent foregrounds and systematic errors faced by the differ- 
ent experiments. The points in Figure^lhave been obtained 
over a frequency range of 26-150 GHz, and by ground-based 
interferometers and balloon-borne scanned total-power tele- 
scopes. They are all from different regions of the sky, the cali- 
brations are all independent and based on different absolute 
calibration sources, and for the two low- frequency experi- 
ments foreground sources have been subtracted in different 
ways, yet the agreement of the power spectra is striking. 

A detailed comparison between the experiments is dif- 
ficult to do from the data points alone because of the corre- 
lated errors between points for each experiment. To make a 
meaningful comparison it is necessary to fit the underlying 
power spectrum to each data set, taking into account the 
correlations, and to compare the parameters that describe 
the power spectrum. In principle many parameterisations of 
the power spectrum would suffice for this comparison, but in 
practice it is obviously sensible to use the standard adiabatic 
cold dark matter power spectrum models, and fits to these 
models for the VSA and other experiments are considered 
in detail in Paper IV. 



Evolution of Stars and Galaxies", under contract HPMT- 
CT-2000-00132. K. Maisinger acknowledges support from 
an EU Marie Curie Fellowship. A. Slosar acknowledges the 
support of St. Johns College, Cambridge. We thank Profes- 
sor Jasper Wall for assistance and advice throughout the 
project. 



7 CONCLUSIONS 

We have derived the power spectrum of the CMB 
anisotropics from the first year's VSA observations, made 
using its compact array configuration. We measure the flat 
band power in 10 weakly-correlated bins of width = 83 
between i ~ 150 and i ~ 900. The results are subject 
to a calibration uncertainty of ±7 percent in power, with 
negligible beam uncertainty. The contribution to the power 
spectrum from diffuse Galactic emission and residual radio 
sources is also negligible. Our results are in excellent agree- 
ment with other recent measurements as regards the ampli- 
tude and position of two peaks in the power spectrum; power 
is also detected out to the resolution limit of the experiment 
at £ ~ 900. 

Band powers, correlation matrices and window func- 
tions are available from 

http : //mrao . cam. ac .uk/telescopes/vsa/results .html. 
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Figure 5. Combined CMB power spectrum from the three mosaiced VSA fields. The error-bars represent Icr limits; the two sets of data 
points correspond to alternative interleaved binnings of the data. 
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Figure 6. A comparison of the VSA data (filled eircles) with results from the BOOMERANG (open squares), MAXIMA (open cireles) 
and DASI (open triangles) experiments. Two sets of error bars are plotted for each data set; the smaller of the two indicate only random 
errors, whilst the larger bars indicate the amount by which the inner points could move due to absolute calibration and beam uncertainty. 
In each case the error bars indicate Icr limits. 
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